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ABSTRACT 

Abell 3667 is the archetype of a merging cluster with radio relics. The NW radio relic is the brightest 
cluster relic or halo known, and is believed to be due to a strong merger shock. We have observed 
the NW relic for ~ 40 ksec of net XMM time. We observe a global decline of temperature across the 
relic from 6 to 1 keV, similar to the Suzaku results. Our new observations reveal a sharp change of 
both temperature and surface brightness near the position of the relic. The increased X-ray emission 
on the relic can be equivalently well described by either a thermal or nonthermal spectral model. The 
parameters of the thermal model are consistent with a Mach number M. ~ 2 shock and a shock speed 
of ^1200 km s _1 . The energy content of the relativistic particles in the radio relic can be explained 
if they are (re)-accelerated by the shock with an efficiency of ~0.2%.. Comparing the limit on the 
inverse Compton X-ray emission with the measured radio synchrotron emission, we set a lower limit 
to the magnetic field in the relic of 3 [iG. If the emission from the relic is non-thermal, this lower 
limit is in fact the required magnetic field. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Abell 3667) — inter- 
galactic medium — radio continuum: general — shock waves — X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Many clusters of galaxies appear to be forming at 
the present time through massive cluster mergers. Ma- 
jor cluster mergers are the most energetic events which 
have occurred since the Big Bang, involving total en- 
ergies of ~10 64 ergs. Merger shocks driven into the 
intracluster gas are the primary heating mechanism of 
the gas in massive clusters. Chandra and XMM-Newton 
X-ray observations have provided beautiful images and 
spectra of merger hydrodynamical effe cts, including cold 
fronts and merger bow shocks (e.g., iMarkevitch et all 
l2000l [20021 ). However, the number of clusters with well- 
observed merger shocks is limited. 

Enigmatic extended cluster radio sources with very 
steep spectra and no clear optical counte rparts have been 
known for ov er 30 years l|Willsonl 119761 for a review see 
lFerettill2008t ). These sources occur in a small fraction of 
rich galaxy clusters. They have such steep radio spectra 
that they can be detected only at lower frequencies (<5 
GHz). Fairly symmetric sources which are projected on 
the cluster center are of ten referred to as "radio halos" 
(e.g., the Coma cluster; iDeiss et aL| [l997). while similar 
elongated sources usually located on the cluster periph- 
ery are called "relics." The extremely steep radio spectra 
flatten somewhat at low frequencies, which suggest that 
the relativistic electrons have undergone significant losses 
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due to synchrotron and inverse Compton (IC) emission. 

In every case, such diffuse cluster radio sources have 
been found in irregular clusters which are apparently un- 
dergoing mergers. This suggests that the radio emitting 
electrons are accelerated or re-accelerated by shocks or 
turbulence associated with these cluster mergers. One 
possible theoretical picture is that the radio halos are ac- 
celerated by turbulence following the passage of merger 
shocks, while the relics are the dire ct result of merger 
shock acceleration (e.g., lFerettill2008l ). 

The same relativistic electrons which produce the radio 
synchrotron radiation will produce hard X-ray (HXR) 
emission by inverse Compton scattering of Cosmic Mi- 
crowave Background (CMB) photons. However, the de- 
tection of this emission has been difficult and remains 
controversial (e.g., comp are iFusco-Femiano et all [2004; 
iRossetti fc Molendi|[200l . 

In many ways, Abell 3667 is the ideal site to study 
mergers and radio relics. It is a very bright X-ray 
cluster at a low redshift (z = 0.0552). The ROSAT 
and ASCA observations showe d that it is a spectacula r 
merger with shock heated gas (jMarkevitch et aflll999h . 
The galaxy spatia l and redshift dist ributions are also 
strongly bimodal (jOwers et all 120091 ) . A3667 was one 
of the first clusters i n which a merger cold front was 
found with Chandra (jVikhlinin et al . 2001b. a|), and the 
term was coined based on the A3667 observations. The 
very sharp density and temperature gradients in this 
cold front have provided very strong limits on the role 
of transport processes (like heat conduction) in clus- 
ters. Subsequent Chandra and XMM/Newto n observa- 
tions (|Mazzotta et all 120021 : [Briel et all 120041 ) have pro- 
vided spectacularly detailed images and information on 
the dynamics of the merger. However, most of the Chan- 
dra and XMM observations study the interior regions of 
the cluster, not the outer regions near the radio relics. 

Abell 3667 contains a pair of curved cluster radio relics 
(|Rottgering et al1 ll997). Their location on either side of 
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Fig. 1. — XMM-Newton mosaic image of A3667 in the 0.5-2 keV energy band. The gain in spatial coverage of the cluster due to the 
addition of our new northwest observation is evident. The image has been background subtracted and exposure corrected. 



the cluster center at 1.8Mpc distance (29') and sharp, in- 
wardly curved outer edges are exactly what is expected 
for merger shocks and shock particle acceleration. Mod- 
els which reproduce the optical and inner X-ray proper- 
ties of A3667 predict shocks at or near these locations 
(jRoettiger et all 119991: iRicker fc Sarazinll2001l) . The ra- 
dio spectra steepen wi th distance from the outer edge 
(|R5ttgering et al.lll!)!)7; as expected if the electrons are 
accelerated there, and the higher energy electrons lose 
energy due to synchrotron and IC emission as they are 
advected away from the shock. 

The northwest radio relic in Abell 3667 is the brightest 
(highest flux) cluster radio relic or hal o source which is 
know n, with a flux at 20 cm of 3.7 Jy (j Johnston- H ollitt 
2004). It is also one of the largest relics, with a total ex- 
tent of 33' or 2.1 Mpc. Since the electrons which produce 
the IC HXR emission are basically the same ones which 



emit radio synchrotron, for a given magnetic field the 
HXR flux should be nearly proportional to the radio flux. 
Thus, Abell 3667 might be expected to be the brightest 
nonthermal (NT) cluster hard X-ray source. Since this 
relic is at a large projected radius from the cluster cen- 
ter, if anything the magnetic field should be lower than 
in other objects, implying an even larger HXR flux. 

We have a long Suzaku observation of the NW radio 
relic region of Abell 3667 with Suzaku (Nakazaw a et al.l 
2009), as well as two observations of more central re- 
gions of the cluster. While hard X-ray emission is 
seen c oming from the clust er center with HXD intru- 
ment (iNakaza wa et al.l I2009T ) . it could be thermal emis- 
sion, while there is no clear Suzaku detection of inverse 
Compton hard X-ray emission from the NW radio relic, 
given the systematic uncertainties in the non- X-ray back- 
ground and thermal emission from the cluster. Our 
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Fig. 2. — Wavelet reconstruction of the 0.5-2 keV X-ray mosaic of A3667 from XMM. There is a sharp discontinuity in the X-ray surface 
brightness to the northwest; The blue ellipses show the location, name and size of three main galaxy components (KMM5, KMM2 and 
KMM4) of A3667 from the analysis of Owers et al. (2009). The black cross marks the center of KMM2 group. The white contours show 
the wavelet reconstruction of the SUMMS 843 MHz radio image in the region of radio relics. The contours are drawn at the 1, 3, 9 and 20 
mjy/beam levels. 
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Suzaku observation also provides a CCD X-ray image of 
a portion of the radio relic region with the X-ray Imaging 
Spectrometer (XIS). However, this covers only a part of 
the region of the radio relic and possible merger shock. 

Here, we present the results of a new XMM-Newton 
observation of Abell 3667 which is centered on the NW 
radio relic. We assume a concordance cosmology with 
H = 71 km s" 1 Mpc" 1 , Q M = 0.27, and il x = 0.73. 
At the redshift of Abell 3667 of z = 0.05525, the angular 
diameter distance is <1a = 212.6 Mpc, and V corresponds 
to 61.8 kpc. For a mean cluster temperature of 7.2 keV, 
the virial radius should be ri 8 o ~ 2.33 Mpc = 38' (?). 
Unless otherwise specified, we provide confidence inter- 
vals at the 68% level. 

2. OBSERVATION AND DATA ANALYSIS 

There have been seven previous XMM-Newton obser- 
vations of Abell 3667. Six of these were discussed in 
iBriel et ail ()2004f ) , where these observations were merged 
to make a mosaic of the center of the cluster. There 
was a subsequent observation (OBSID 0206850101) of 
the center of the cluster, which we analyzed and added 
to the mosaic for completeness. All these observations, 
however, cover only the central part of the cluster. In or- 
der to study the characteristics of the X-ray emission at 
the location of radio relic, we proposed for an additional 
XMM-Newton observation, with the pointing direction 
selected to ensure a robust background subtraction can 
be performed by leaving enough area beyond the loca- 
tion of relic. The observational ID of this pointing is 
0553180101, and it was performed during XMM-Newton 
orbit 1620. In this paper, we will concentrate on the 
results obtained only from this observation. The total 
scheduled time of observation was 55 ksec. We used 
XMMSAS (Watson et al. 2001) version 7.3 for the stan- 
dard data reduction. After producing the event files, 
we have applied a strict light curve cleaning in order to 
ensure that we retained only the periods with the lowest 
and the most stable background. The resulting net expo- 
sures after the cleaning are 37, 47 and 44 ksec for the pn, 
MOS1, and MOS2 detectors, respectively. The medium 
filter has been used for all instruments, for consistency 
with the previous observations of the cluster. Extended 
full frame mode has been used for the pn to reduce the 
effect of out-of-time events. In producing the mosaic, we 
have checked whether any of the ccds exhibit unusually 
high background, identifying none. 

3. X-RAY IMAGE ANALYSIS 

We us ed the quadruple background subtraction tech- 
nique of iFinoguenov et all ((2007) to produce the mosaic 
image of A3667. Due to high level of cluster emission, 
several iterations have been performed in order to locate 
and mask out the bright parts of the cluster. We used 
both pn and MOS data in producing the mosaic image. 
Figure Q] shows the mosaic image made with the new ob- 
servation and the previous data from XMM in the 0.5-2 
keV band. It has been corrected for background, expo- 
sure, vignetting, and out-of-time events, and excluded 
the 1.4-1.6 keV energy band due to the strong back- 
ground Al line. Compar i son to previous XMM X-ray 
images (jBriel et al.l 120041 : iNakazawa et al.1 l2009f ) shows 
the significant gain in spatial coverage achieved with the 
new observation. The image in Figure [2] shows a cor- 



responding wav e let re construction using the method of 
I Vikhlinin et al.l (|1998h and applying the image recon- 
struction technique of Finoguenov et al. (2009) to remove 
the wings of point sources. 

To study the spatial correlation of the X-ray emis- 
sion to the radio relic, we used the SUMMS radio data. 
SUMSS is a deep radio survey of the entire sky south of 
declination —30 degrees, made using the Molonglo Ob- 
servatory Synthesis Telescope, operati ng at 843 MHz an d 
recording right-circular polarization (|Bock et al.l Il999h . 
The beam size of these radio images is 45". We used 
the data which are a part of SUMSS release 0.2. To 
characterize the radio relics, we apply the wavelet recon- 
struction of radio image, remove point-like features and 
display the contours at the 1, 3, 9 and 20 mJy/beam lev- 
els in Figure [2] The wavelet reconstruction of the XMM 
image of Abell 3667 shows an interesting coincidence be- 
tween the extent of the significant X-ray emission in the 
0.5-2 keV detected on the wavelet scales upto 8' and the 
outer edge of the radio relic. The detection stops 11' 
from the end of the data, compared to 2' on the oppo- 
site edge of the cluster. As discussed in § [IJ this is the 
expected location of the cluster merger shock, if the relic 
was produced by shock acceleration or re-acceleration of 
relativistic electrons. 

Our previous Suzaku XIS image of this region did not 
show this edge (jNakazawa e t al. 2009) , which we believe 
to be a result of contamination by point sources which 
were not individually resolved due to the large Suzaku 
XIS point spread function (PSF). Many of these point 
sources are seen in Figures [T] & [5] In addition there is a 
steep gradient in the cluster emission to the southeast, 
associated with the major subcluster, which is smeared 
by the large Suzaku XIS PSF and covers the relic. To test 
the effects of these features and the Suzaku XIS PSF on 
the Suzaku image, and the consistency of the Suzaku and 
XMM data, we blurred the XMM image using a Gaussian 
with a sigma of 1', which is a crude approximation to 
the Suzaku XIS PSF. We found that this blurred image 
matched the actual Suzaku XIS image fairly well, and 
the surface brightness discontinuity observed with XMM 
at its full resolution is no longer apparent in the smeared 
image. 

A more detailed examination of the XMM image shows 
in addition to the X-ray peak associated with the main 
cluster, another peak, with coordinates of the X-ray cen- 
ter are RA = 20' i 10 m 41 i ;285, Deck = -56°35'45'.'81 
(J2000). The available results of the spectroscopic survey 
of the cluster (|Johnston-Hollitt et aLll2008l lOwers et atl 
l2009h allows us to identify this peak with the major sub - 
cluster in A2667, KMM2, found by (jOwers et all 12009( 1. 
which plays a key role in the dynamics of the system. 
In fact the relics are not symmetric with respect to any 
of the galaxy component but the region between KMM2 
and KMM5, indicating that within the forward shock hy- 
pothesis of their origin, it is the recent passage of KMM2 
through KMM5 that is responsible for their appearence. 
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Fig. 3. — The zones masked out from both the surface brightness 
and spectral analysis. We removed all significant point sources, as 
well as the bulk of the cluster emission and the emission from the 
two northwest galaxy groups. The dashed box shows the region 
for surface brightness extraction presented in Fig|4] The black 
contours show the wavelet reconstruction of the SUMMS 843 MHz 
radio image in the region of radio relics. The contours are drawn 
at the 1, 3, 9 and 20 mjy/beam levels. 

In order to avoid biased results caused by the possibly 
different temperature gas associated with this subcluster, 
we excluded the region from further analysis of both the 
cluster and relic emission, together with all the point and 
compact extended sources detected in the image. We 
used the wavelet images on spatial scales of 4, 8, 16, 32, 
and 64 arcsec to mask out such regions. Figure |3] shows 
all of the regions masked out from the surface brightness 
and spectral analysis. 

In Figure 2J we show the X-ray surface brightness ra- 
dial profile around the radio relic. We used the 1.6-4.0 
keV band to extract the profile and use only pn data. The 
low-energy cut off has been selected to avoid any possible 
influen ce of the strong lo w-energy emission seen towards 
A3667 (|Briel et al.ll2004ft . which is likely of Galactic ori- 
gin. The upper energy is selected to maximize the signal- 
to-noise, as at higher energies the background dominates. 
The background in the selected energy band is dominated 
by particle background which is not vignetted, and thus 
is roughly constant in surface brightness. We see a sharp 
drop in the surface brightness profile at the position of 
the outer edge of the relic, indicating that there is an as- 
sociation of this X-ray edge with the radio relic. A small 
peak in the radio surface brightness at +10' distance from 
the relic is accosiated with two radio point sources. We 
associate a corresponding low-significance enhancement 
in X-rays is two sources below the detection threshold. 
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Fig. 4. — The surface brightness profile of A3667 in the 1.6- 
4.0 keV band extracted along the direction from the cluster center 
(negative distance to the relic) to the northwest. The radial units 
are arcminutcs. The location of the outer edge of the relic is at 
0. Point sources were excised. The dashed line shows the surface 
brightness profile of radio emission, extracted in the same fashion 
from the 843 MHz SUMMS data. One sees a sharp drop in the 
surface brightness at the position of the relic and a much flatter 
behaviour both at small and at large X-axis values. 

In Figure [SJ we show the combined pn+MOS XMM- 
Newton image in the 1.6-4 keV band for the region im- 
mediately around the radio relic. There is a clear asso- 
ciation of the surface brightness discontinuity with the 
location of the outer edge of radio relic. This image 
was not corrected for vignetting to show that the result 
is not affected by the details of background subtraction 
(vignetting correction in combination with undersubtrac- 
tion of the background can lead to enhancements at the 
edge of the detector) . The image is smoothed with a 32" 
Gaussian after the removal of point sources. The X-ray 
image appears to follow the radio relic in some detail. 
In particular, the X-ray surface brightness discontinuity 
traces the outer, sharp edge of the radio relic, including 
the indentation at the west end of the relic. To quantify 
the trends seen in Figure [SJ we proceed with the spectral 
analysis. 

4. SPECTRAL ANALYSIS 

In order to investigate the origin of the break in the 
surface brightness, we extracted the X-ray spectra from 
several regions around the radio relic. The regions were 
designed to encompass the zone directly in front of the 
relic, much of the relic itself, as well as other regions trac- 
ing the possible azimuthal variations, and cluster radial 
temperature profile. These spectral extraction regions 
are shown in Figure [6l We used the spectral range of 
0.4-14 keV in the spectral analysis. We started by fitting 
the most distant zone (from the cluster center) covered 
by XMM observation. The energies above 7 keV where 
clearly dominated by the residual particle background, 
which we fit with a background (not convolved with the 
effective area of the telescope) power law model. The 
residual emission was found to consists of a kT ~ 0.3 
keV thermal emission we associate with the Galactic fore- 
ground emisson and a significant 1 keV component, con- 
sistent with the Suzaku XIS results, which we therefore 
associate with the A3667 cluster emission. The outmost 
zone is also the largest zone and the constraints on the 
background components are the best there, so we fixed 
the slope of the particle background and the tempera- 
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Fig. 5. — XMM-Newton image (pn+MOS) of the region around the NW radio relic in A3667 in the 1.6-4 keV energy band. The X-ray 
point sources have been removed. The image was smoothed with a 32" Gaussian. The white contours show the radio emission without the 
wavelet reconstruction applied. Dark zones indicate the regions of intense X-ray emission. The image shows that the complex shape of the 
radio emission is followed by the X-rays. 

ture of the soft emission to the values found there, while 
the normalization of both these components was left free 
to vary. We found this background/foreground model to 
be a satisfactory fit, and the reported additional thermal 
and non-thermal components we associate with the clus- 
ter are clearly in excess of it. Since no spatial variation 
in the particle background spectral shape has been re- 
ported [except for the case of bright MOS chips (Kunz 
& Snowden 2008), which we associated with electronic 
noise and not with particle background itself], we are 
certain that the reported results correspond to a detec- 
tion of a change in the intrinsic properties of the cluster « *" ® § 
emission. 
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Fig. 6. — The spectral extraction regions for the NW radio relic. 
The numbers correspond to the Zone IDs in Tables [T] & [2] 

The results for the spectral fits for the cluster emis- 
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TABLE 1 

Thermal Model Spectral Fits to the X-ray Emission near the Relic 





kT 


Norm. a 






Area 


I x (2-5 keV) 




Zone ID 


fkeV) 


do -5 ) 


v 2 

A r 




(arcmin 2 ) 


(10 — ergs/s/cm^ /arcmin^) 


Comment 


1 


1.0 ±0.2 


1.0 ±0.2 


0.90 


297 


60.4 


0.9 ±0.6 


background - main 


2 


4.0 ±0.8 


17.5 ± 1.7 


0.90 


111 


23.3 


8.9 ± 1.0 


relic back-1 


3 


4.7+g'l 

— U.o 


15 ± 1.5 


0.82 


75 


21.8 


13.2 ± 1.8 


relic back-2 




k 1+4.1 
S - i -1.6 


Q 1 A- 1 1 




149 


102 1 


a q -j- n ^ 
*±.o m u.o 


v r 
relic iront 


5 


5 4+ 3 ' 3 


6.0 ±0.6 


1.33 


54 


28.0 


8.0 ± 1.3 


relic back-3 


6 


3.8 ±0.8 


24 ± 1.4 


0.96 


131 


80.8 


11.4 ± 1.2 


cluster near relic 


7 


6.0 ± 1.3 


26 ± 1.7 


0.89 


91 


63.2 


19.2 ± 1.8 


cluster near relic 


8 


4.7 ± 1.2 


24 ± 1.7 


1.10 


100 


82.4 


20.4 ±2.1 


cluster near relic 


9 


1.0 ±0.1 


4.1 ± 0.5 


1.19 


19 


56.9 


2.2 ±0.8 


zone in front of relic 


10 


1.9 ±0.6 


12 ±2 


1.03 


441 


765.2 


0.9 ±0.3 


zone in front of relic 


11 


1.7±0.8 


6.9 ± 1.7 


1.06 


137 


259.2 


0.6 ± 1.0 


zone in front of relic 


12 


2.5 ±0.5 


37 ±2.2 


1.05 


290 


503.6 


4.3 ±0.7 


eastern relic 


13 


1.8 ±0.6 


15 ±2 


0.91 


245 


476.6 


1.7 ±0.5 


zone west of relic 


14 


1.0 ±0.3 


0.12 ±0.05 


0.92 


136 


340.3 


2.2 ± 1.2 


background side 


15 


1.6 ±0.4 


4.8 ± 2 


1.38 


19 


75.7 


6.0 ± 1.5 


zone behind the relic 



a Normalization of the APEC thermal spectrum, which is given by {10 _14 /[47r(l + z) 2 ^^]} J n e riH dV, where z is the 
redshift, (La is the angular diameter distance, n e is the electron density, ng is the ionized hydrogen density, and V is 
the volume of the region. 



TABLE 2 

Power-Law Spectral Fits to the Relic X-ray Emission 



Zone ID 


r 


Norm. a 


X r 


N d .o. f . 


comment 


2 


1.6 ± 0.1 


15 ±4 


0.91 


Ill 


relic back-1 


3 


1.8 ±0.1 


38 ±2 


0.79 


75 


relic back-2 


1 


1.4 ±0.2 


17 ±3 


1.09 


149 


relic front 


5 


1.9 ± 0.2 


16 ± 2 


1.28 


51 


relic back-3 



a The normalization of the power-law is the photon flux at 1 
keV in units of 10 — 6 photons cm -2 s _1 keV -1 . 



sion components are given in Tables Q] & [2j The first 
table gives the result of purely thermal spectral fits us- 
ing APEC plasma code (Smith et al. 2001), while the fits 
in Table [5] include nonthermal power-law components. 

We first consider the purely thermal fits to the spec- 
tra. The results are summarized in Table [TJ where the 
columns give the Zone ID from Figure [51 the tempera- 
ture T, the normalization of the APEC thermal model, 
the reduced chi-square Xr, the number of degrees of free- 
dom in the fit Nd.o.f., the projected solid angle of the 
region excluding the area masked to avoid point sources, 
the average X-ray surface brightness of the region Sx in 
the 1.6-4.0 keV band, and a comment on the reason for 
choosing this region. The element abundance has been 
set to 0.2 solar (typical for cluster outskirts), as no strong 
lines have been seen in the emission to constrain it. 




Fig. 7. — The reconstructed surface brightness in the 2-5 keV 
band using the results of the spectral fits, and directly fitting for the 
spectral shape of the background without assuming a flat spatial 
distribution. The grey scale shows the flux level in units of ergs 
s" 1 cur 2 arcmin -2 . We note that an enhancement in surface 
brightness is located on the relic, as seen by the difference in regions 
on either side of the relic. The larger circle shown is centered 
on the cluster center. The smaller circle is centered on the main 
subcluster. 

As one moves inward in radius towards the relic, there 
is initially a slow rise in the temperature, from 1 to 2 keV, 
followed by an abrupt increase in the temperature from 
kT ~ 1.9 keV in front of the relic to ~ 5 keV on the relic. 
Figure [7] shows the average X-ray surface brightnesses of 
each of the regions plotted along with contours showing 
the radio relic. We confirm the dent-like structure of the 
X-ray emission, seen in FigJSJ Figures [4] & [7] both show 
that there is also a jump in the X-ray surface brightness 
inward of the outer edge of the radio relic. The value of 
the jump is higher in Figure [3 compared to Figure [U 
which we attribute to the presence of the dent in the X- 
ray emission, which was separated though the selection 
of spectral extraction zones, but not in the extraction of 
the surface brightness profile. 

The combination of the jump in the temperature and 
the jump in the surface brightness (which is determined 
mainly by the gas density) are consistent with a location 
of the shock at or near the front of the relic. The param- 
eters of the shock are discussed in the next section. In 
any case, whatever its origin, the hardening of the spec- 
trum at the front of the relic is significant at better than 
the 95% confidence level. 

Alternatively, the jump in the X-ray surface brightness 
and spectral hardening might also be attributed to the 
non-thermal emission associated with the relic. To ad- 
dress this possibility, we also fit the spectra in the regions 
corresponding to the projected position of the NW radio 
relic with a non-thermal power-law spectrum. The re- 
sults of these fits are summarized in Table [2j where the 
columns give the Zone ID from Figure O the power-law 
photon spectral index T, the normalization of the power- 
law which is the flux at 1 keV, the reduced chi-square Xr, 
the number of degrees of freedom in the fit Nd.o.f., and 
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the same comment on the location of the region as given 
in Table [T] Note that the photon spectral index values 
are generally lower than might have been expected from 
the steep radio spectrum of the radio relic. The average 
radio spectral index bet ween 85 MHz and 2.3 G Hz for 
the whole relic is a » 1.1 (Rott gering et al1ll997| h which 
corresponds to V w 2.1. However, the fits in Tabled] 
are for pure power-law emission. Almost certainly, there 
is also some thermal emission from the intracluster gas, 
which may affect the fitted spectral index of the power- 
law. The limited statistics in the observed X-ray spectra 
of these regions do not allow for more complicated ther- 
mal plus power-law models. Also, for reasonable values 
of the magnetic field, the X-ray emission in the relatively 
soft XMM band would come from lower energy relativis- 
tic electrons than produce most of the observed radio 
emission. The curvature in the observed radio spectrum 
indicates that the electron energy spectrum flattens at 
lower energies. A steep radio spectrum at high frequen- 
cies is expected as a consequence of the radiative losses 
by the relativistic electrons. 

Unfortunately, as a result of the limited statistics in the 
spectra of these faint regions in the outer parts of this 
cluster, it isn't possible at this time to determine whether 
thermal or non-thermal models provide a better fit to the 
data. If one considers zone 4 as an example, the thermal 
model gave a very slightly better fit, with a reduced chi- 
square of \l = 1-06, as compared to the power-law model 
which had \r = 1-09- The x 2_ test indicates that both 
models are acceptable. The difference becomes even less 
significant if one combines the fits to zones 2-5. 

5. DISCUSSION 

With the new XMM data on the regions around the 
northwest radio relic in Abell 3667, we have found a dis- 
continuity in the X-ray surface brightness at the outer 
edge of the radio relic. This increase in the X-ray sur- 
face brightness inside of the relic compared to outside 
also corresponds to an abrupt hardening in the X-ray 
spectrum. If the X-ray emission from the relic is thermal 
emission from intracluster gas, then this feature corre- 
sponds to a jump in the density, temperature, pressure, 
and specific entropy of the gas. The most likely origin of 
such a jump would be a merger shock in the cluster. 

One can estimate the shock Mach number M. from the 
compression behind the shock. From the X-ray surface 
brightness on either side of the surface brightness discon- 
tinuity, we estimate that the pre-shock and post-shock 
electron densities are n e i = (6.81 ± 0.55) x 10~ 5 cm -3 
and n ea = (1.32 ±0.08) x 10~ 4 cm -3 , respectively. This 
gives a shock compression of C = n e , 2 /n e i = 1.94 ±0.19, 
where the subscripts 1 and 2 refer to pre-shock and post- 
shock gas. For a 7 = 5/3 shock, the shock jump condi- 
tions give 



C AM 2 4' y ' 

This gives a shock Mach number of Al = 1.68 ± 0.16. 
The effects of projection, the finite width of the regions 
used to determine the surface brightness, the XMM PSF, 
and the likely curvature of the shock front mean that this 
is probably an underestimate. 

One can also estimate the Mach number from the tem- 
perature jump at the shock front. Using the values for 
zones 10 and 4 for the pre-shock and post-shock gas, we 



have a temperature jump of T2/T1 = 2.68 ± 1.19. For a 
7 = 5/3 shock, the jump conditions give 

T2 _ 5M 4 + UM 2 - 3 

7i ~ 16M 2 ' ( ' 

Using the observed temperature jump and propagating 
the errors gives M = 2.43 ± 0.77. Because of the larger 
errors in the temperature, the uncertainty in this value 
is large. 

If one takes the weighted mean of the Mach number 
from the shock compression and shock temperature front, 
one finds Ai = 1.71 ± 0.16. Using this value, the shock 
compression is C = 1.97 ± 0.19. The pre-shock sound 
speed is c s \ = 710 ± 110 km s _1 which gives a shock 
speed of v s = M.c s \ = 1210±220 km s _1 . Merger shocks 
of similar and slightly higher strength has been detected 
in the bullet cluster and A520 (Markevitch et al. 2002; 
2005). However, the location of the shocks in these clus- 
ters is in the core region, while in A3667 the location 
of the shock is in the low-density region. Perhaps, the 
only other cluster where the shock has been found in the 
low-density region is A754 (Krivonos et al. 2003; Henry 
et al. 2004). 

Given the very low gas densities in this region roughly 
2.2 Mpc from the center of the cluster, and the high shock 
speed, it is possible that the post-shock gas has not had 
time to come into electron-ion equipartition by Coulomb 
collisions, or into collisional ionization equilibrium. If 
the electrons are not strongly heated by the shock, the 
time scale for the electron and ion te mperatures to come 
into equipartition is a pproximately (jFox fc Loebl 119971 : 
iWong fc SarazTnll27)09h 

«■"•"*«• 8 (ifef (if^)"W'- 

where T e is the electron temperature, n p is the proton 
density, and In A is the Coulomb logarithm. Because the 
possible merger shock is weak (has a low Mach num- 
ber), the electrons will be heated significantly by adi- 
abatic compression, even if there is no shock electron 
heating. For the adapted value for the shock compres- 
sion, adiabatic heating gives (T e2 /T e i) s» 1.57, while the 
full shock heating including adiabatic compression gives 
(T 2 /Ti) ks 1.72. If the adapted post-shock density and 
the post-shock electron temperature assuming only adi- 
abatic heating are used in equation (|3]), the approximate 
time to reach equipartition is t cq fs 3 x 10 s yr. The 
speed of the post-shock material relative to the shock 
is t>2 = v s /C « 610 km s~ x . Thus, the thickness of 
the region with a lowered electron temperature would 
be d eq ~ V2t e q ~ 0.2 Mpc, corresponding to an angular 
scale of 9 cq ss 3cos</> arcmin. Here, where <f> is the angle 
between the central shock normal and the plane of the 
sky. This is comparable to the width of our spectral zone 
4. So, this would argue that nonequipartition could be 
quite important. On the other hand, the apparent post- 
shock temperature is actually a bit higher than expected 
assuming full shock heating and the shock compression 
derived from the jump in X-ray surface brightness. This 
may indicate that the electrons are heated effectively at 
the shock. 

At the observed post-shock temperature, one would ex- 
pect the gas to achieve ionization equilibrium after a time 
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iionoq ~3x 10 12 /n e sec (e.g.. iFuiita et al. I [20081. For 
the post-shock electron density, this would give tioncq ~ 
7 x 10 s yr, which is even longer than the timescale to 
reach equipartition. However, the effects of nonequilib- 
rium ionization may be difficult to discern without much 
better spectra than we currently have. 

One argument favoring a thermal shock model for the 
jump in brightness and hardness of the X-rays at the 
outer edge of the relic is that a similar jump is seen in 
nearby regions where the radio emission from the relic 
is not bright. For example, the zones which are closest 
to the apparent shock position but in regions with weak 
radio emission are zone 2 for the shocked and zone 13 
for the preshock. The density jump between these two 
zones corresponds to a Mach number of M. « 2.2. The 
temperature jump is basically identical to that measured 
on the radio relic, so the estimate of Mach number from 
the temperature jump would not change. Thus, the jump 
in X-ray properties is similar in regions with bright or 
faint radio emission, which is consistent with thermal 
emission from a shock, but perhaps not with IC emission 
from the relic. 

Could this merger shock be accelerating or re- 
accelerating the relativistic electrons in the NW radio 
relic? We start by estimating the energy dissipated in 
the shock. As a simple estimate of this, we take the 
change in the kinetic energy flux across the shock. This 
is given by 



AR 



1 



KE 



1 - 



1 

c 2 



(4) 



where p\ is the pre-shock mass density in the gas. Us- 
ing our derived values for the shock properties gives 
Ai^KE ~ 9 x 10~ 5 ergem -2 s _1 . The width of the relic 
from northeast to southwest is roughly 26.'3 or 1.63 Mpc. 
Taking the area of the shock perpendicular to the flow 
as a circle with this diameter gives a perpendicular area 
of 2.09 Mpc 2 . With this size, the total rate of conversion 
of shock kinetic energy is 

^ « 1.8 xlO 45 ergs- 1 . (5) 
at 

We now compare the rate dissipation in the shock to 
the energy required to power the radio source. Given the 
steep radio spectrum of the relic, it appears that the rel- 
ativistic electrons are losing energy due to radiation. We 
first consider the losses due to radio synchrotron emis- 
sion. Assuming the relic is in steady-state, with the shock 
accelerating new relativistic electrons and these particles 
losing energy by radiation, the total rate of energy de- 
position in the electrons is given by the radio luminosity. 
Unfortunately, there are few recent radio observations of 
the relic, and some of the existing older values disagree, 
suggesting that some flux in some observations may be 
missing due to the large extent and low surface bright- 
ness of the relic, or that some of the flu xes include other 
sourc es. The flux at 1.4 GHz is 3.7 Jy ([Johnston-H ollitt 
l2004f) . and the overal l spectral index from 8 5.5 MHz to 
2.4 GHz is a » -1.1 (|Rottgering et al.lll997| ). The spec- 
tral index appears to steepen to a « —1.9 from 1.4 GHz 
to 2.4 GHz, and to flatten below 1.4 GHz. Thus, we 
approximate the radio spectrum as 



-0.9 



-1.9 



. 1.4GHz/ 



v < 1.4 GHz 
v > 1.4 GHz 



(6) 



Then, the total radio luminosity of the relic is L ra dio ~ 
3.8 x 10 42 erg s _1 . The IC luminosity of the relic is given 
by Lie = ^radio(3.6^G/S) 2 , so the total nonthermal 
luminosity is L NT « 3.8 x 10 42 [(3.6 nG/B) 2 + 1] erg s" 1 . 
If this energy is provided by the shock acceleration of 
electrons, then the efficiency of shock acceleration is 



d.E c 
dt 
<lt: he 
dt 



dE K E 
dt 



0.0021 



3.6 fiG 

B 



1 



(7) 



Below, we show that the magnetic field is the relic is at 
least B > 3 pG, so that the correction factor in equa- 
tion ([7]) is between one and 2.5. Note that this acceler- 
ation efficiency is about one order of magnitude smaller 
than the values of a few percent usually inferred from 
the radio emission by Gala ctic supernova remnants (e.g., 
iRosswog fc B rugger] 120071 ). This might be due to the 
relatively low Mach number of this shock compared to 
these in supernova remnants. Alternatively, it may be 
that part of the X-ray discontinuity at outer edge of the 
radio relic is due to IC emission, and the shock is actually 
weaker than our estimate. 

First order Fermi acceleration gives relativistic elec- 
trons with a power-law spectrum n{E) dE oc E~ p dE, 
where the power-law index is p = (C + 2)/(C — 1) and 
C is the shock compression. The spectral index for ra- 
dio emission near the shock should be a = — (p — l)/2, 
while the integrated radio spectrum of the relic and spec- 
trum of hard X-ray inverse Compton emission should be 
one power steeper. Using the value of the compression 
determined above, this would give a w —1.55 at the 
shock, and a « —2.55 for the integrated spectrum of 
the radio relic. The observed radio spectra are flatter, 
with a value of a ss —0.7 near the outer edge of the 
relic, and —1.9 < a < —1.1 for the inte grated spectrum, 
depending on the observed frequency ([Rott gorin g et all 
119971 Uohnston-Hollittl 120041 ). This might indicate that 
we have underestimated the shock compression, or that 
the simplest shock acceleration model for the relic does 
not apply. 

The fact that the radio spectrum steepens with 
projected distance from the outer edge of the relic 
(|R5ttgering et al.l Il997f ) , and that the relic fades 
away to the southwest rather than ending abruptly 
(jRottgering et al.lll997tlJohnston-Hollittll2004() might be 
explained if the relativistic electrons are accelerated at 
the shock, and then undergo radiative losses as they are 
advected away from the shock. If the radio spectrum 
steepens at a frequency Uf,, the radiative age for elec- 
trons which produce the break is (jvan der Laan fc Perolal 
1969) 



trad « 1-3 X 10 8 (- 
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-1/2 



B 
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(8) 



The speed of the post-shock material relative to the shock 
is V2 — v s /C « 610 km s _1 . Thus, the relativistic elec- 
trons will have moved a distance c? ra d away from the 
shock, where 



^rad ' 



i 0.082 



1.4 GHz 
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The corresponding angular distance is 



Mpc . 
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rad ' 
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1.4 GHz 



-1/2 



3.6 ^G 
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-i -l 



-3/2 
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COS ( 



(10) 



where <f> is the angle between the central shock normal 
and the plane of the sky. This is roughly consistent with 
the observed thickness of the layer at the front outer edge 
of the radio relic where the spectrum b etween 20 and 13 
cm is observed to steepen dramatically (Rottgcrin g et al] 
1997). On the other hand, the full width of the relic is 
about 8' from NW to SE at 1.4 GHz. Since this corre- 
sponds to several electron radiation loss lengths, this may 
indicate that there is re-acceleration of particles within 
the relic, perhaps by turbulence produced by the pas- 
sage of the merger shock. Alternatively, far from the 
NW edge of the relic, we may be seeing radio emission 
from relativistic electrons which have been recently ac- 
celerated and which are located at the front or back edge 
of a convex shock region . 

This argument neglected the diffusion of the relativis- 
tic electrons. If Bohm diffusion ap plies, the diffusion 
coefficient is (Bcrczins kv et al. 



ion ap r. 

Mm, 



D(E) w 2 x 10 



29 



E \ 1/3 
lGeV J 



B 

TJIg 



-1/3 



2 — 1 

cm s 



(ii) 

The electrons which produce the radio emission at 20 
cm have E « 10(5/3 /iG) -1 / 2 GeV, which leads to 
D(E) w 3 x 10 29 (B/3nG)-^ 2 cm 2 s" 1 . During the 
radiative loss time t la d, these electrons would diffuse a 



distance of dam ~ [D(E)t Ta d 



1/2 



0.01 Mpc, which is 



smaller than the advection distance. We also note that 
the sharp outer edge of the radio relic suggests that dif- 
fusion is not very important on the scale of the size of 
the relic. 

Alternatively, the increase in the X-ray brightness and 
hardening of the spectrum at the outer edge of the relic 
might be the result of inverse Compton emission, rather 
than thermal emission from the shocked gas. As note in 
§ SI the signal-to-noise of the X-ray spectra do not al- 
low these two models to be distinguished. If the jump in 
X-ray emission and hardness is due to inverse Compton 
emission, one would still need to explain why the relic 
is located where it is. That is, how are the relativistic 
electrons in the relic generated? It might be that the par- 
ticles are due to first order Fermi re (acceleration) from 
a merger shock; in this case, there would still need to be 
a shock at the front of the relic, but the X-ray emission 
in this region might be predominantly nonthermal. It is 
also possible that the electrons are accelerated by some 
other mechanism. 

Taking the nonthermal fits as an upper limit on the 
true flux of inverse Compton emission coming from the 
relic in regions 4 and 5, we can derive a lower limit 
on the average magn e tic fie ld strength B as shown by 
lHarris <§Fll omanishm1 (|1974l ). Considering monochro- 
matic fluxes Fr(i>r), Fx(vx) at frequencies vr, vx, for 



a power-law distribution of electrons, the expression for 
the magnetic field is 



B = C(p)(l + z)^l^+V ( |k 

„ \ (P-1)/(P+1) 

vx ) 



2/(P+l) 



(12) 



The value of the proportionality constant C(p) can be 
found from the ratio of the sync hrotron flux (lLongair 
1994 eqn. 18.49) to the IC flux (jRvbicki fc Light man 
19791 eqn. 7.31). For region 4 we use the photon in- 
dex found for region 5, since this is more consistent with 
the spectral index of the radio relic. The lower pho- 
ton index found in region 4 may be the result of coin- 
cident thermal emission, and therefore our estimate of 
B would be low. Assuming the X-ray flux in these re- 
gions is nonthermal with a photon index in the range 
1-9 ^ r < 2, we find magnetic field strengths in the range 
2.3 /iG < B < 4.5 fj,G. Since it is unknown what frac- 
tion of the emission is thermal, B may in fact be larger. 
Therefore we estimate a lower limit for the magnetic field 
at the location of the radio relic of B > 3 /iG. This limit is 
just consistent with Faraday rotation measure estimates 
from the the observation o f background radio gala xies in 
the zone close to the relic (|Johnston-Hollitdl2004| ). 

6. CONCLUSIONS 

We have made a new X-ray observation of the north- 
west radio relic region in the Abell 3667 cluster in order 
to gain insight on both thermal and non-thermal phe- 
nomena that should be associated with it. We detect a 
jump in both the surface brightness and the hardness of 
the X-ray emission at the outer edge of the relic. The 
X-ray spectra of these regions is consistent with the ex- 
tra emission on the relic being either thermal emission 
from a merger shock or inverse Compton emission from 
relativistic electrons in the radio relic. There are some 
arguments that favor the shock model, which implies a 
shock with a Mach number of M. w 2 and a shock speed 
of ~T200 km s _1 . The energy content of the relativis- 
tic particles in the radio relic can be explained if they 
are (re)-accelerated by the shock with an efficiency of 
~0.2%. The merger shock model and shock acceleration 
of the relativistic electrons can explain, at least roughly, 
the radio surface brightness and spectral distributions in 
the relic. 

Alternatively, the jump in the brightness and hardness 
of the X-rays at the outer edge of the relic might be due 
to inverse Compton emission from the relic. Assigning 
all the emission from the relic to the IC component, we 
derive a lower limit on the magnetic field in the relic of 
>3 /iG. 

These observations have yielded the first robust char- 
acterization of the X-ray properties near the northwest 
radio relic and possible merger shock in Abell 3667. Un- 
fortunately, the data are not sufficiently deep to deter- 
mine if the relative roles the thermal shock emission and 
nonthermal IC emission. Recently, we were approved 
for a very long observation of this regions as an XMM- 
Newton Large Program. These data should allow us to 
determine the efficiency of particle acceleration and the 
relative roles of thermal and nonthermal processes in this 
merger shock. 
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